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Abstract 

Beta-integers ("/3 -integers") are those numbers which are the counterparts of integers when 
real numbers are expressed in irrational basis (5 > 1. In quasicrystalline studies j3 -integers 
supersede the "crystallographic" ordinary integers. When the number )3 is a Parry number, 
the corresponding )3 -integers realize only a finite number of distances between consecutive 
elements and somewhat appear like ordinary integers, mainly in an asymptotic sense. In this 
letter we make precise this asymptotic behavior by proving four theorems concerning Parry 
j3 -integers. 

1 Introduction: j5 -integers versus integers in quasicrystals 

Aperiodicity of quasicrystals fl\ implies the absence of space-group symmetry to which integers 
are inherent 0. On the other hand, experimentally observed quasicrystals show self-similarity (e.g. 
in their diffraction pattern). The involved self-similarity factors are quadratic Pisot-Vijayaraghavan 
(PV) units. We recall that an algebraic integer j8 > 1 is called a PV number if all its Galois conju- 
gates have modulus less than one. Mostly observed factors are the golden mean x = and its 
square = , associated to decagonal and icosahedral quasicrystals respectively. The other 
ones are 5 = 1 + \/2 (octagonal symmetry) and 8 = 2 + \/3 (dodecagonal symmetry). Each of 
these numbers, here generically denoted by j8, determines a discrete set of real numbers, Z^g = 
{bn I « G Z} fS\. The set of j8 -integers bn aims to play the role that integers play in crystallog- 
raphy. In the same spirit, j8 -lattices Y are based on j3 -integers, like lattices are based on integers: 
r = Ya^i I'pGi, with (e,) a base of R^. These point sets in M'^ are eligible frames in which one could 
think of the properties of quasiperiodic point-sets and tilings, thus generalizing the notion of lattice 
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in periodic cases ||4l. As a matter of fact, it has become like a paradigm that geometrical supports of 
quasi-crystalline structures should be Delaunay sets obtained through "Cut-and-Projection" from 
higher-dimensional lattices. Now, it can be proved Q that most of such Cut and Project sets are 
subsets of suitably rescaled j3 -lattices. 

j8 -integers and their j3 -lattice extensions have interesting arithmetic and diffractive properties 
which generalize to some extent the additive, multiplicative and diffractive properties of integers 
lHIH H. Throughout these features, a question arises to which the content of this paper will 
partially answer. In view of the "quasi"-periodic distribution of the ft„'s on the real line, it is 
natural to investigate the extent to which they differ from ordinary integers, according to the nature 
of j8. If j8 is endowed with specific algebraic properties, one will see that the corresponding set of 
j8 -integers is a mathematically "well-controlled" perturbation of Z in an asymptotic sense 

bn ~ CR{n + a{n)) (1) 

where « i-^ is a bounded sequence and the scaling coefficient is in Q(j8). As a matter of 
fact, for a large family of numbers j8, namely the Parry family, j3 -lattices asymptotically appear 
as lattices. Parry numbers give rise to j3 -integers which realize only a finite number of distances 
between consecutive elements and so appear as the most comparable to ordinary integers. 

We give in the present letter a set of rigorous results concerning the asymptotic behavior of jS- 
integers when j3 is a generic Parry number. More precisely, we will illustrate the similarity between 
sets N and = {bn | « G N} for j8 being a Parry number by proving two properties: 

bn 

1 . We will show that cr = lim„^oo — exists and we will provide a simple formula for cn . 

n 

2. For j8 being moreover a PV number with mutually distinct roots of its Parry polynomial, we 
will prove that {bn — cp n)„eN is a bounded sequence. 

Let us mention that both of the previous asymptotic characteristics are known for j8 being 
a quadratic unit. The following proposition providing explicit formulae for jS -integers for a quadratic 
unit /3 comes from 171 . 

Proposition 1.1, If fi is a simple quadratic Parry unit, then 

= {^^^ = + ^ 1^ + ^ { ^ } > « e n} , Where c, = ^1±^ . 
If P is a non-simple quadratic Parry unit, then 

where {x} designates the fractional part of a nonnegative real number x. 
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In Section 2, we provide a brief but self-contained description of the concept of j8 -integers, 
with emphasis on jS -integers realizing only a finite number of distances between neighbors, i.e., 
/3-integers associated with Parry numbers j3. The distances are then coded by letters, and, thanks 
to the self-similarity of Zjg, the obtained infinite words up are known to be fixed points of sub- 
stitutions. As the associated substitution matrices are primitive, the Perron-Frobenius theorem 
together with some further applications from matrix theory will lead to the results in Section 3 on 
the asymptotic behavior of j3 -integers for Parry numbers j3. 



2 All we need to know on j8 -integers 
2.1 p -representation and j3 -expansion 

Let j3 > 1 be a real number and let xhe a non-negative real number. Any convergent series of the 

form X = Y!i=-o^Xil5' XkXk-\ ■ ■ -xq^x-i ■ ■ ■ , where x,- G N = {0, 1,2, . .. } and x^ ^ 0, is called 
a p -representation of x. Let j8 be a positive integer, then if we admit only {0, 1, . . . , j8 — 1} as 
the set of coefficients and if we avoid the suffix (j3 — 1)®, where ft) signifies an infinite repetition, 
there exists a unique j3 -representation for every x, called the standard j3 -representation. For j3 = 10 
(resp. P = 2), it is the usual decimal (resp. binary) representation. 

Even if /3 is not an integer, every positive number x has at least one j3 -representation. This 
representation can be obtained by the following greedy algorithm: 



1. Find keZ such that j3* < x < j3*+i and put Xk := 
X — {x} denotes the lower integer part of x. 

2. For i < k, putx, := Li3'";+iJ and := {j3r,+i}. 



X 



pkj andrfe:= |^ }>, where [x\ 



This representation is called the P -expansion of x and the coefficients of the jS -expansion 
clearly satisfy: x^ G {1, • • • , [P\} and xi G {0, . . . , [j3J} for all / < k. We use the notation {x)p 
for the j3 -expansion of x. For j3 being an integer, the j3 -expansion coincides with the standard 
P -representation. 

The greedy algorithm implies that among /3 -representations, the j3 -expansion is the largest 
according to the radix order. Actually the latter corresponds to the ordering of real numbers: for all 
x,y e [0,-|-oo), the inequaUty x <y holds if and only if (x)j3 is smaller than (y)jg according to the 
radix order. 

Example 2.1. Let jS = T = ^^y^. The golden mean T is the larger root of the polynomial x^ —x—l. 
Applying the greedy algorithm, we get for instance the following 1 -expansions: 
(^) =0.1, (3^) =100., (5^) =i.(oooir---. 
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2.2 Renyi expansion of unity 

The j8-expansion of numbers from the interval [0, 1) can be obtained through the map Tp : [0, 1] 
[0,1) defined by 

Tp{x)={l5x}. (2) 
It is easy to verify that for every x G [0,1), it holds (x)^ = 0»x_ix_2''' if and only if x_,- = 
[pTp ^'^^^(x)J . By extending Formula ^ to x = I, one gets the so-called Renyi expansion of unity 
m in base j8: 

dp{\)=tit2t3--- , where f,- := [jSr^^^ (1)J . (3) 

Every number j3 > 1 is characterized by its Renyi expansion of unity. Note that ti = [j8J > 1. 
Parry in |9| has moreover shown that the Renyi expansion of unity enables us to decide whether 
a given j8 -representation of x is its j3 -expansion or not. For this purpose, we define the infinite 
Renyi expansion of unity (it is the largest infinite j3 -representation of 1 with respect to the radix 
order) 

/ if dp{\) is infinite, 

' 1 (?lf2 -I))'" if dp{\)=ty-t„, with t^^O. ^' 

Proposition 2.2 (Parry condition). Let dp{l) be the infinite Renyi expansion of unity in base j8. Let 

Y!i=-oo^iP' be a ^-representation of a non-negative number x. Then Lf=-oo-^!j8' is the ^-expansion 
ofx if and only if 

XiXi-\ ■■■ -< d*p{l) for all i < k, (5) 
where -< means smaller in the lexicographical order. 

Example 2.3. For j3 = T = , the Renyi expansion of unity is <iT(l) = 11. Then, d*{\) = (10)®, 
and, according to the Parry condition, any sequence of coefficients in {0, 1} which does not end 
with (10)'" and which does not contain the block 11 is the z-expansion of a non-negative real 
number 



2.3 Parry numbers 

This paper is concerned by real numbers j3 > 1 having an eventually periodic Renyi expansion 
of unity dp{\), i.e., dp{\) = t\ . . .tm{tm+\ ■ ■ - tm+p)^ for some m,p G N, called Parry numbers. 
For every Parry number j8, it is easy to recover, from the eventual periodicity of dp{\), a monic 
polynomial with integer coefficients having j8 as a root, i.e., j8 is an algebraic integer. However, this 
so-called Parry polynomial is not necessarily the minimal polynomial of /3. It was proven in lUOll 
that every Pisot number is a Parry number. 

2.4 Definition and properties of /3 -integers 

Nonnegative numbers x with vanishing j8 -fractional part are called nonnegative ^-integers, for- 
mally, := {■'c > I (x)^ = XkXk-\ •••xo»}. The set of 15-integers is then defined by := 
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(— Z^) UZp. As already mentioned, the radix order on j3-expansions corresponds to the natu- 
ral order of non-negative real numbers. Consequently, there exists a strictly increasing sequence 
(^n)r=o such that 

^0 = and {bn\n£N} = Z+. (6) 

When j8 is an integer, = Z and so the distance between the neighboring elements of Zp is 
always 1. The situation changes dramatically if j3 N. In this case, the number of different 
distances between the neighboring elements of Zp is at least 2 and the set Zp keeps only partially 
the resemblance to Z: 

1. Zp has no accumulation points. 

2. Zp is relatively dense, i.e., the distances between consecutive elements of Zp are bounded. 

3. Z|3 is self-similar, thus pZp cZp. 

4. Zp is not invariant under translation. 

5. Zp forms a Meyer set if j8 is a Pisot number, i.e., Zp — Zp C Zp + F for a finite set C M 
(proved in ||3l). 

Thurston ifTTTl has shown that distances occurring between neighbors ofZp form the set {A^^ | G N}, 
where 

^-t-L^ for^GN. (7) 
(=1 P 

It is evident that the set {A^. | /c G N} is finite if and only ii dp (l) is eventually periodic. 
2.5 Infinite words and substitutions associated with /3 -integers 

If the number of distances between neighbors in Zp is finite, one can associate with every distance 
a letter. Thus, we obtain an infinite word up coding Z'^ as illustrated in Figure [U In order to define 

= 010010100100101... 

10 

• • • • • 

It ^2^j 

Figure 1: First elements of Z+ (non-negative T-integers) and of the associated infinite word u^. 

Up properly, let us introduce some definitions from combinatorics on words. An alphabet is 
a finite set of symbols, called letters. A concatenation of letters is a word. The set s^* of all finite 
words (including the empty word £) provided with the operation of concatenation is a free monoid. 
The length of a word w = wqWiWi • • • w„_i is denoted by |w| = n and \w\a denotes the number of 
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letters ae^inw. We will deal also with infinite words u = uqu\U2 ■■■ & A finite word w 
is called a factor of the word u (finite or infinite) if there exist a finite word w^^^ and a word w^^^ 
(finite or infinite) such that u = w^^^ww^^K The word w is a prefix of u if w*^'^ = £. A substitution 
over ^* is a morphism <p : ^* s^* such that there exists a letter a ^ and a non-empty 
word w G s^* satisfying 9(0) = aw and <p is non-erasing, i.e., 9(^7) / £ for all b ^ s^. Since any 
morphism satisfies <p(vw) = (p{v)(p{w) for all v,w G jz/*, the substitution is uniquely determined 
by images of letters. Instead of classical (p{a) = w, we sometimes write a^w.K substitution can 
be naturally extended to infinite words u G s^'^^ by the prescription (p{u) = (p{uo)(p{ui)(p{u2) . . . 
An infinite word u is said to be a fixed point of the substitution <p if it fulfills u = (p{u). It is 
obvious that a substitution <p has at least one fixed point, namely lim„^oo (p"{a)- Morphisms on j^* 
equipped with the operation of composition form a monoid. Substitutions do not form a monoid as 
the identity is not a substitution. Nevertheless, if <p is a substitution, 9" is also a substitution. 

The following prescription associates with every substitution <p defined over the alphabet ^ = 
{ai , a2 , . . . , } a non-negative integer dxd matrix, called the substitution matrix M^, 

{M^).. = \(p{ai)\aj, (8) 

As an immediate consequence of the definition, it holds for any word w that 

= (|<P Wla, , |<P W|a„ . . . , |<P WlaJ. (9) 

The substitution matrix of the composition of substitutions 9, i/a obeys the formula M,po =^v/^<P- 

Example 2.4. The Fibonacci substitution (p defined on {0,1} by (p{0) = 0\, 9(1) = 0, has the 

following substitution matrix 

Let us finally recall a notion from matrix theory which proves useful in the study of substitu- 
tions. A matrix M is called primitive if there exists ^ G N such that all entries of M*^ are positive. 
There exists a powerful theorem treating primitive matrices. 

Theorem 2.5 (Perron-Frobenius). Let M be a dxd primitive matrix. Then: 

1. the matrix M has a positive eigenvalue X which is strictly greater than the modulus of any 
other eigenvalue, 

2. the eigenvalue X is algebraically simple, 

3. to this eigenvalue corresponds a positive eigenvector (i.e., with positive entries only), while 
no other eigenvalue has a positive eigenvector. 

The eigenvalue X from the above theorem is called the Perron-Frobenius eigenvalue of M. If 
a primitive matrix M is the substitution matrix of a substitution 9 over jz/ = {ai , . . . according 
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to the result of Queffelec |[T2l . the left eigenvector (pi , p2, . . . , Pd) of A normalized by ^f^j p, = 1 
is equal to the vector of letter frequencies in any fixed point u of (p, i.e., 

(pi, p2,...,pd) = (p(ai), p{a2),...,piad)), (10) 

where p(ai) = lim„^oo '- and uS"' denotes the prefix of u of length n. 

n 

2.6 Parry numbers and the infinite words 

From the formula for distances we know that the number of distances in Z^g is finite if and only 
if the Renyi expansion of unity d^{\) is eventually periodic, i.e., if jS is a Parry number. 

• If dp{\) is finite, i.e., dp{\)= t\t2 ■ ■ ■ t,,,, t^^O, p is said to be a simple Parry number, and the 
set of distances is {Aq, Ai, . . . ,A,„_i}, where all of the listed elements are mutually distinct. 

• If (jfj3 (1) is eventually periodic, but not finite, j8 is a non-simple Parry number. Choose p,m^ 
N to be minimal such that dp{\) = ?i?2 • • -^mCfm+i • • •^m+p)'*'. then the set of all mutually 
distinct distances is {Aq, Ai, . . . , Am+p-i}. 

Let us precisely define the infinite word u^ = uqU\U2 ■ ■ ■ associated with for a Parry number 
j3. Let {A(), . . . , Aj_i} be the set of distances between neighboring j3 -integers and let {b„)'^^Q be as 
defined in ©, then 

def 

Un = i if bn+l-bn = Ai. (11) 

2.7 Canonical substitutions for Parry numbers 

Fabre in lfT3l has associated with Parry numbers canonical substitutions in the following way. 

Let jS be a simple Parry number, i.e., dp{l) = t\t2 ■■■tm, for m E N. Then the corresponding 
canonical substitution (p is defined over the alphabet {0, 1, ... ,m — 1} by 

(p(0) =0''l, (p(l) =0'22, ... ,(p(m-2) =0'"-'(m-l), (p(m- 1) = 0'"'. (12) 

Similarly, let j8 be a non-simple Parry number, i.e., dp{l) = t\t2 ■ ■ -tmitm+i ■ ■ ■tm+p)"' ■ The associ- 
ated canonical substitution (p is defined over the alphabet {0, 1 , . . . , m + p — 1} by 

9(0) =0''l, (p(l) =0'^2, ... ,(p(m- 1) =0''"m, ... 

... ,(p{m + p-2) =0''"+"-'{m + p- 1), (p{m + p-\) = 0'"^"m. (13) 

Each of these substitutions has a unique fixed point lim„^oo <p"(0). Moreover, this fixed point turns 
out to be equal to up. It is readily seen that in both cases, the substitution matrix Mip is primitive. 
The characteristic polynomial of the substitution matrix M(p coincides with the Parry polynomial 
of the number j3. Hence, the Perron-Frobenius theorem implies that the Parry number jS is a simple 
root of its Parry polynomial p{x) and therefore p'{p) / 0. 
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3 Asymptotic behavior of jS -integers for Parry numbers 



In order to derive some information about asymptotic properties of j8 -integers, let us recall the 
essential relation given in [.13.1 between a j8 -integer bn and its coding by a prefix of the associated 
infinite word M^g . 

Proposition 3.1. Let up be the infinite word associated with a Parry number (5 and let (p be the 
associated substitution of jS, then for €V€ry hyi G it holds that (J^n)^ — ^k—v • • 'CIiciq^ if dud 
only if (p^-^{0"t-' )... (p(0"' )0"° is a prefix of up of length n. 

Since every prefix of up codes a j8 -integer bn, Proposition 13.11 provides us with the following 
corollary. 

Corollary 3.2. Let w be a prefix of up, then there exist ^ G N and ao, ay,. .. ,ak-i G N such that 

w = (p*^-i(0"*-')...<p(0'")0''°, 
where at-i .. .aiao* is the P -expansion of a p -integer. 

Let us denote Ui := |(p'(0)|, then (f/,)Jlo ^ canonical numeration system associated with the 
Parry number /3 (defined and called j8- numeration system in lITOl and studied in fl3l ). For details 
on numeration systems consult fT4l. Proposition 13. II implies that the greedy representation of an 
integer n in this system is given by 

k-i k-i 
n=Y, aiUi if bn = Y^ a,j8'. 

i=0 (=0 

Applying dS), the sequence {Ui)Y^Q may be expressed employing the substitution matrix M of cp in 
the following way 

[/, = (1,0,...,0)M'(1,1,...,1)^ (14) 

where T is for matrix transposition. Let j8 be a simple Parry number, then the Renyi expansion 
of unity is of the form dp{l) =tit2...tm and j8 is the largest root of the Parry polynomial p{x) = 
x"' — {tix"'^^ +f2x'"^^ H hfm-i-^ + fm)- Let US recall that p{x) may be reducible. 

Our first aim is to provide a simple formula for constant cp such that b„ « cpn. For any root 

7 of p{x), it is easy to verify that (/"^^ Y"^^, • • • , 7, 1) is a left eigenvector of the substitution 
matrix M associated with 7. On the other hand, according to (flOl ). the unique left eigenvector 
(po , Pi Pm- 1 ) of M with Pi = 1 is such that p,- is the frequency of letter / in up . Combining 
the two previous facts, we obtain for frequencies the following formula 



Let (Ao, Ai, . . . ,Am-i)^ be the right eigenvector of M associated with jS such that Aq = 1, then 
it is easy to verify that A, is the distance between consecutive j8 -integers which is coded by letter 
/ in the infinite word up (see the formula (O for distances). For our purposes, the following easily 
derivable formula for distances will be useful 

A- = i8' - L 0/3'"^ /G{0,l,...,m-1}. (16) 

7=1 

Theorem 3.3. Let p{x) be the Parry polynomial of a simple Parry number j3. Then 

b„ P-l 

Proof. Let us denote by u the prefix of up of length n, then 

bn = |m|oAo + |m|iAi H h \u\m-l^m-l- 

Since frequencies of letters exist, lim„^oo ^ exists and obeys the following formula 

b„ 

hm — =poAo + piAiH hp,«-iA,„_i. 

Applying (ITSl) and (fT6l ). we obtain 



lim„_^ = 1 (l-i^-i-XiS'-L^^iOiS'"^')) 
^ - I-r/ tj{m - j)^-'-^ 



□ 

Corollary 3.4. Let j8 = jSi, j82, • • • , jS,,, Z7e mutually distinct roots of the Parry polynomial p{x) of 
a simple Parry number jS. Then 



^ m 

lim 



i=2 

Proof p{x) = Y\7=i - Pi), p'{x) = iti nr=i,,-^.(^ - A), thus p'{p)= nr=2(i8 - ao- □ 

Remark 3.5. If p{x) is an irreducible polynomial, then jS is an algebraic integer of order m and 
p2,--- ,Pm are algebraic conjugates of j5, and hence mutually distinct. 
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We now study the asymptotic behavior of {bn — cpn)ne'H- We know already that the limit 
lim„^oo ^ exists. Hence it is enough to consider the limit of the subsequence 

lim — = lim = lim — . 

Under the assumption that all roots of p{x) are mutually different, we will find a useful expression 
for Un- Since M is diagonalizable, there exists an invertible matrix P such that PMP^^ is diagonal 
with {PMP^^).. = pi, / G {1, . . . ,m}. Using ([HI), we may write 

Un = ih 0,...,0)p-l(PMP-l)>(l,l,...,l)^ (17) 

It follows from the Perron-Frobenius theorem that j8 > |j8, |, hence, the formula (T7\ leads to the 
following expression 

- = lim^ = (l, 0,...,0)P-1P,,P(1,1,...,1)^ (18) 

where Pg, is the orthogonal projection on the line given by ei = (1,0, . . . ,0)^, i.e., {Pe,)ii = 
1, {Pe,)ij = otherwise. We now examine the difference b„ — cpn. Let {b,^)^ = at-i ...ao», 
thus b„ = L/=o '^iP' ^ — L?Jo Employing ([TtI) and (fTSl) . we obtain 

-b„ - n = £ «Ki , 0, . . . , o)p-' (/3;p,. - (PMp-i) > (1 , 1 , • • • , 1)^ = 

= (1,0,...,0)P-1ZP(1,1,...,1)^ (19) 

where Z is a diagonal matrix with Zn =0, Zjj = —zj for j G {2, . . . ,m}, and zj = <^iP'j- Since 
the coefficients of j8-expansion satisfy a,- G {0,..., [jSJ} and since for PV numbers j8, it holds 
IjSyl < 1 for j = 2,3,... ,m, we have 

Remark 3.6. Suppose that the Parry polynomial p(x) of a Parry number p is reducible, say p{x) = 
q{x) ■ r{x), where q{x) is the minimal polynomial of fi, and r{x) is a polynomial of degree at least 1. 
Then the product of the roots of r{x) is an integer and therefore either all roots of r{x) lie on the 
unit circle or at least one among the roots ofr{x) is in modulus larger than 1. It implies that the set 
ofzj is bounded for all j if and only if fi is a Pisot number and its Parry polynomial is the minimal 
polynomial of (5. 

According to Remark [T6l and as P does not depend on n, we have shown the following theorem. 
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Theorem 3.7. Let (5 be a simple Parry number If j3 is moreover a Pisot number and the Parry 
polynomial of p is its minimal polynomial, then {b„ — cp is a bounded sequence. 

Now f is a matrix of the Vandermonde's type given by 

/^, = A'"-,then Pil,l,...,lf = (^,^^,...,mY. ill) 



In order to have for all n e N an exphcit formula for -^b„ — n, it remains to determine ( 1 , 0, . . . , 0)P^ ^ , 

i.e., the first row of p-K Since P^ = ^P'"^K where {P"''^)ij = detP{j, 1) and P(j, 1) is obtained 
from P by deleting the j'-th row and the 1-st column, applying Vandermonde's result yields 



.p-u _ n,<Ki,kMi^i-h) i-iy-' i-iy-' .... 



Notice that since p{x) does not have multiple roots, p'{^j) 7^ 0. 
It follows from formulae OH), (HI]), and that 

1, 



In consequence, using the estimate (l20l) . we may deduce an upper bound on \bn — cp n\ 

"'1 1 
\K-cpn\< 2cpPl^ , ■ (23) 



Example 3.8. Let us illustrate the previous results on the case of the simplest simple Parry number 
- the golden mean j8 = T = Renyi expansion of unity is d-z{\) = 11 and p{x) = x^ —x — I. The 
substitution matrix for the Fibonacci substitution has been given in Example \2.4\ Consequently, 
iUn)nef'i satisfies Un = fnfor all « G N, where {fn)nen is the Fibonacci sequence given by 

fn+\ = fn+fn-l, /o = 1, /l = 2. 

Applying Theorem \3. 3\ we get 

_ T-1 , _ p'{z) _ 2t-1 _ T^ + 1 

which is in correspondence with Proposition li.il 

Let us denote the second root of p{x) (the Galois conjugate of T) by t', x' = ^—j-^ = If 
{bn)p = ak-i ■ --aiao; then at G {0, 1} and 

, I 1 , / \ / 1 + T \ l-x . 
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Since \ < 1, the sequence \b„ — «|nGN bounded and we may easily determine an upper bound 
(taking into account that t' < Oj 

thus, comparing the upper and lower bound, we deduce 

\bn-c^n\ < 



which is again in correspondence with Proposition \l.l\ where we have replaced the fractional part 
{f+^} ^^^^ ^ in order to get an upper bound on |Z7„— Cf«|- 

3.1 Non-simple Parry numbers /3 

Let j8 be a non-simple Parry number, then the Renyi expansion of unity is of the form <ij3 ( 1 ) = 
hh ■ ■ -tmitm+i ■ ■ -tm+p)^ with m,p chosen to be minimal and j3 is the largest root of the Parry poly- 
nomial p{x) = {xP — I) (x*" — fix™^' tm) — t,n+\xP^^ — tm+p-\x — t^+p- Let US recall 

that p{x) may be reducible. 

Similarly to the simple case, our first goal is to derive a simple formula for constant Cj^ such 
that b„ Can. For any root 7 of the Parry polynomial p{x), 



y-\YP-i),...,(f-i),f-\...,y, 1) 



m components p components 

is a left eigenvector of the substitution matrix M associated with 7. On the other hand, according 
to (fTOl) . the unique left eigenvector (po,pi, . . . ,Pm+p-i) of j8 such that 1^"!^^^ p,- = 1 satisfies that 
Pi is the frequency of letter / in up. Combining the two previous facts, we obtain for frequencies 
the following formula 

p. = ^ = ^'(^-1) (24) 

ES'i8'(i8/'-l) + LfrJi8' )8"'()8P-1)' 

where 

CJi = p'"-'^-'{pP - 1) for < / < m - 1, 
Oj = j3'«+P-i-' for m<i<m + p-l. 

Let (Ao, Ai, . . . , A,„+p_i)^ be the right eigenvector of M associated with j3 such that Aq = 1. 
Then it is easy to verify that A, is the distance between consecutive j3 -integers which is coded by 
letter / in the infinite word up (see the formula for distances in ([Til). Similarly as for simple Parry 
numbers, also for non-simple Parry numbers the following formula for distances holds and will be 
useful 

Ai = p'-'£tjl5'-J, /G{0,l,...,m + p-l}. (25) 
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Theorem 3.9. Let p{x) be the Parry polynomial of the non-simple Parry number j8. Then 

y b„ 15-1 f 

b„ 

Proof. Similarly as for simple Parry numbers, lim„^co — exists and we have 



lim — =poAo + piAiH \- Pm+p-i\n+p-i- 

Applying (l24l) and (|25] |. we obtain lim„^oo — = „ f „ — ^ — r (A +8), where 

ff J fe 1 n ^ j8'"(j8/'-l) 

m-l i m+p-l i 

A=l^r-'-\^'-W-"Ltjr^) and B= £ ^^P-'-^^' -l^tj^^)- 

i=0 7=1 i=m j=l 

It is then straightforward to prove that 

A={^P- 1) (mi8'"-i - tj{m - jW'-^-') , 

A + B = p'{P). 

□ 

Remark 3.10. If we consider the infinite Renyi expansion of unity d*^{\) instead of the "classical" 

Renyi expansion of unity dp{l), we have in the simple Parry case dp{l) = (fi . ■■tm-\{tm ~ I))"- 
Thus the length I of preperiod is and the length L of period is m. In the non-simple Parry case, 
we have d*^{\) = dp{\) = titj- ■ - tmitm+i ■ ■ - tm+p)"^- Hence the length I of preperiod is m and the 
length L of period is p. With this notation, the formulae for cpfrom Theorems \3. 3\ and \3. 9\ mav be 
rewritten for both simple and non-simple Parry numbers in a unique way as 

Corollary 3.11. Let jS = j8i,j82, • • • ,Pm+p be mutually different roots of the Parry polynomial p{x) 
of a non-simple Parry number fi. Then, 

CO = lim ^ = , , n (j8 - A). 

Proof. Analogous as in Corollary 13.41 □ 
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Remark 3.12. If p{x) is an irreducible polynomial, then (5 is an algebraic integer of order m + p 
and j82, . . . ,Pm+p <^f^ algebraic conjugates of p, and hence mutually different. 

Let us now investigate the asymptotic behavior of [b^ — cp n)neN- As we know akeady that the 
Umit hm„^c„ ^ exists, we may rewrite it in terms of the subsequence {U„) 

Um — = Um = Um — . 

Under the assumption that all roots of p{x) are mutually distinct, we will express {]„ in an easier 
form. Since M is diagonalizable, there exists an invertible P such that PMP^^ is diagonal with 
[PMP^^).. = jii, i£ {\,...,m + p}. Using (O, we may write 

Un = {l, 0,...,0)p-l(PMP-l)'>(l,l,...,l)^ (26) 

It follows from the Perron-Frobenius theorem that j8 > |j8, |, hence, the formula (l26l ) leads to the 
following expression 

- = lim^ = (l, 0,...,0)P-1P„P(1,1,...,1)^ (27) 

Now, let us turn our attention to the difference b„ — cpn. Let {b„)p = ak-\ ■ --aQ*, thus bn = 
Ef=o ^i^^ ^^'^ " = I^to ^i^i- Employing (l26l) and (l27l) . we obtain 

1 k-\ 

-bn-n=Y^ai{\,0,...,0)p-\P\P,,-{PMP-'y)P{l,l,...,lf = 

= (1,0,...,0)P-1ZP(1,1,...,1)^ (28) 

where Z is a diagonal matrix with Zn =0, Zyy = —zj for 7 G {2, . . . ,m + p}, and zj = L/=o '^iPj- 
Since the coefficients of j3 -expansion satisfy a,- G {0, . . . , [j3J } and since for PV numbers p, it holds 
IjSyl < 1 for j = 2,3,... ,m + p, we have 

\zA<tm\<-^y (29) 

According to Remark [3^ and since P does not depend on n, we have shown the following theorem. 

Theorem 3.13. Let (5 be a non-simple Parry number If /3 is moreover a Pisot number and the 
Parry polynomial of [5 is its minimal polynomial, then {b„ — cp n)ne^ is a bounded sequence. 
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The explicit form of the matrix P reads 



1) 
1) 



3m 



(Pi- 



^2 



\ A 



Hence, 



m-lrop 
'm+p Kr'm+p 



pfia. 



om-2(oP 
Hm+p KHm+p 



1) 



(A" 



m+p 



13P-I 
■'m+p 



->m+p 



1 \ 



1 / 



,1^ 



rm-Vp \rm+p 



(30) 



In order to have for all n G N an explicit formula for -^bn — n, it remains to determine ( 1 , 0, . . . , 0)P ' , 

i.e., the first row of P^. By contrast to the simple Parry case, the matrix P is not in the Vander- 
monde's form. However, we notice that its determinant is equal to a Vandermonde determinant 
through a simple addition of columns. More precisely, we start with the addition of the last column 
to the m-th column, the last but one column to the (m — l)-st column and so forth. It is readily 
seen that this procedure leads after m steps to a Vandermonde matrix of order m + p with the same 
determinant as P. 

So detP = Y[i<k{Pi — Pk) ■ The expression of {P^^)ij, j e {I, . . . ,m + p},h then given by 



n 



(A-A^ 



1; 



i<k. i.k^j 

n,<.(A-A^ 



(31) 



Notice that since p{x) does not have multiple roots, p'{Pj) / 0. 
We obtain applying expressions ( |28] ). ( [30l ). and (|3T| ) 



1 



m+p 



In consequence, we may deduce an upper bound on \bn —cpn 

1 



n < 



^''^^S(1-|A;I)^ Ip'iPj) 



(32) 



Example 3.14. Let us illustrate the previous results for the simplest non-simple Parry number j8 
with Renyi expansion of unity dp{\) = 21® and p{x) = x^ — 3x + \, i.e., j8 = = The sub- 

stitution matrix for the associated substitution (p : — > 001, 1 ^ 01 is the square of the Fibonacci 

substitution matrix, i.e., M^p = ( ^ ^ j. Consequently, {Un)neN is just a subsequence of the 

Fibonacci sequence {fn)nen (defined in Example 
Applying Theorem \3.9\ we get 



cp 



1 



A(A-i 



given by Un = fm 
p'm 215-3 



1 



A 



1 ^2' 



15 



which is in correspondence with Proposition li.il 

Let us denote the second root of p{x) by p', j8' 
a,- G {0, 1,2} and 



3-V5 



2 



Let {bn)p = at-i . . .aiao; then 




Since < j8' < 1, the sequence {bn — cp is bounded and we may easily determine an upper 

bound (taking into account that coefficients in (5 -expansions satisfy the Parry condition) 



which is in correspondence with the estimate we get if we replace the fractional part |^| with 1 
in Proposition li.il 

Open questions: 

It would be interesting to study the behaviour of j3 -integers even for non-Pairy numbers j3, 
i.e., when the Renyi expansion of unity of j(3 is not eventually periodic. Among non-Parry 
numbers one may distinguish two cases: 

• If the length of blocks of zero's in cfj3 (1) is bounded, say by a length L, then the j8 -integers 
form a Delone set since the shortest distance between consecutive points is at least and 
the largest distance is 1. 

• If d^{\) contains strings of zero's of unbounded length, then the set of distances (A/t) 
between consecutive jS -integers have as its accumulation point, see Equation (O. It means 
that in this case, the j8 -integers do not form a Delone set. 

The first question to be answered in both cases is whether the hmit bn/n does exist for some 
of non-Parry numbers j8 . 
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